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INTERACTION OF LARGE, HIGH POWER SYSTEMS WITH 
OPERATIONAL ORBIT CHARGED PARTICLE ENVIRONMENTS 

by Carolyn K. Purvis, 'N. John Stevens, 
and ^ Frank D. Borkopee 

Concepts are presently being advanced for space systems 
to be used for such activities as manufacturing, Earth ob- 
servations, scientific exploration, power generation and 
human habitation, in locations ranging from low Earth or- 
bit (900 - 500 km) to geosynchronous orbit and beyond. 

Many of these systems concepts envision large structures 
and high power levels, and consequently higher operating 
voltages than have been used in space to date. The poten- 
tial impact of interactions of space systems with their op- 
erational orbit charged particle environments on the sys- 
tems' performance must be accounted for in the design 
process. A potentially hazardous spacecraft- environment 
interaction is discussed in this paper. This is the inter- 
action of large high voltage systems with low energy (<50 
eV) plasmas which can result in loss of power and/or arc- 
ing. The impact of this class of interactions on system 
operation is most severe at low orbits where the ambient 
plasmas are densest. Results of experimental work and 
predictions of simple analytical models are presented and 
their implications for design of space systems are dis- 
cussed. 

INTRODUCTION 

Studies are presently being conducted for space systems to be used for such di- 
verse activities as manufacturing, scientific exploration, power generation and 
human habitation in locations ranging from low Earth orbit (500 - 400 km) to 
geosynchronous orbit and beyond * • * * . Many concepts for such systems en- 
vision large complex structures and high power levels. Consequently, for 
higher electrical efficiency, higher operating voltages than those used to 
date are being proposed. In designing these new systems, care must be taken 
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to consider the influence of interactions with orbital charged particle environ- 
ments on a system's operation. These systems will be expensive to build and 
launch; their economic feasibility is strongly dependent on optimal design and 
long term reliability. This report presents a preliminary evaluation of the 
effects of charged particle interactions with space systems in Earth orbits 
ranging from shuttle altitudes to geosynchronous altitude (fig. 1). 

The interactions of space systems with their charged particle environments can 
be divided into three general types, with each comprising an interaction of the 
space system with pai’ticles of a characteristic energy level. 

The first of these types is the interaction of high voltage systems with the 
thermal plasmas (temperature <2 eV) of the plasmasphere. These plasmas 
range in density, composition and temperature as a function of altitude. Typical 
random electron current and ram ion current (resulting from spacecraft motion) 
as functions of altitude are shown in figure 2 . This type of interaction is due to 
the voltage differences among various portions of the space system. Spacecraft 
flown to date have had low voltage power systems. Therefore the impact on 
their performance has been negligible. However, for the new generation of high 
power high voltage systems now being proposed, these interactions must be con- 
sidered. 

The primary emphasis of this paper is on the high voltage system interactions. 
The concept of the basic mechanism involved is discussed, preliminary results 
of experiments with a solar array segment are presented, and a simple model 
of a large solar array in low Earth orbit and in geosynchronous orbit is used in 
a preliminary evaluation of the impact of these interactions on the performance 
of a space system. 

The second general type of interaction is the charging of spacecraft in geosyn- 
chronous orbit by kilowatt particles (~1 to ~50 keV)(see fig. 2). During geo- 
magnetic substorm activity, clouds of kilovolt electrons and protons are in- 
jected into the midnight sector of the magnetosphere, and are observed at geo- 
synchronous orbit . These particles are known to charge spacecraft negatively 
with respect to plasma potential to kilovolt levels in eclipse, and to hundreds of 

n 

volts in sunlight . Spacecraft charging and its impact on space systems is the 

g 

subject of an ongoing investigation . There is a substantial literature on space- 
craft charging ’ * * and it will not be discussed further here. 
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The final category of interactions is the radiation damage to spacecraft compo- 
nents caused by hi$i energy particles (~100 keV and greater). This type of in- 
teraction is a known phenomenon and will not be discussed further. 

HIGH VOLTAGE SYSTEM/THERMAL ENVIRONMENT INTERACTIONS 

Basic Interaction Mechanism Concept 

12 

It is well known from probe theory “ that any surface placed in a plasma will 
collect current from the plasma and that this collection is dependent upon the 
surface voltage. Any system immersed in a plasma will come into equilibrium 
with the plasma by acquiring the charges necessary to make the net current be- 
tween the system and the plasma zero. 

In general, the electron flux (and therefore electron current density) in a plasma 
is larger than the ion flux. A conducting surface (simple probe) placed in a 
plasma will thus collect more electrons than ions. This will cause it to take on 
a negative voltage witli respect to the plasma. It will then collect fewer elec- 
trons and more ions. In equilibrium it will have acquired negative potential 
such that the electron and ion currents to it are equal, i. e. it is collecting net 
zero current. Typically this negative "floating'" potential will be of order the 
electron energy. 

In the case of a system with a distributed voltage, a solar array for example, the 
system will float at voltages relative to the plasma so that the positive area col- 
lects electrons and the negative area collects an equal current of ions. Fig- 
ure 3 depicts this concept for a solar array having ;ui operating voltage V 
with a linear distribution along the panel length. The negative portion of the 
array collects ions, and the positive portion electrons. Since, in general, elec- 
tron fluxes are higher than ion fluxes in a plasma, a larger area is required for 
collecting ions. Consequently, in equilibrium, the magnitude of the negative 
voltage (collecting ions) will be larger than that of the positive voltage. 


The net effect of the process of particle collection is that a current (I T 
passes through the plasma effectively creating a parasitic load (Rpj \,sm \) ’ n 
parallel with the electrical load of the array. This results in a loss of useable 


array power. 
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Experimental Results 


To calculate the magnitude of this effect, it is necessary to make an estimate 
of the current collection of an array which has small biased conductor areas 
surrounded by large areas of insulation. This estimate is made with the help of 
data from an experiment using a biased array segment. The blowup in figure 3 
shows the construction of a current lightweight large solar array design. The 
sketch is a cross section through a single cell. The cells are 2x4 cm and 
0.02 cm (8 mils) thick with wraparound contacts exposed only on the edges, and 
are covered by 0. 015 cm (6 mil) fused silica cover glasses. They are mounted 
on a substrate consisting of copper interconnects sandwiched between two layers 
of 0. 001 cm (0. 5 mil) kapton. The cell connectors are welded to the intercon- 
nects through holes in the kapton; there are access holes to the copper inter- 
connects in the back (non cell) side of the array. Thus, while each surface of 
the array is largely covered with insulators, there is some exposed metal on 
each side. 

A segment of an array of this construction has been tested in the Lewis Research 
Center's (LeRC) geomagnetic substorm simulation facility 1- *. The segment, ob- 
tained from the Marshall Space Flight Center, is 14 cm x 20 cm (see fig. 4) and 
has 20 cells mounted in a series-parallel circuit. A plasma was generated by 
ionizing nitrogen gas in the facility's plasma source. The pi'operties of the en- 
vironment ior Uie tests conducted on this segment were: ambient pressure 
~10 _<> torr, plasma density ~10 "Vem'* and electron temperature ~3 eV. The 
electrical circuit of the array was biased to -*] KV in a series of steps by 
means of a power supply outside the tank. The voltage profiles across the array 
segment were obtained by sweeping the facility electrostatic probe across the 
sample at fixed time interval. Thus, transient to steady- state profiles were ob- 
tained. The tests were conducted twice, once with the probe surveying the cell 
side and secondly with the substrate lx’ing surveyed. All test results presented 
here were obtained without sunlight and with the sample nominally at room 
temperature. Prior testing has verified that sunlight and temperature do not 
cause significant changes in the test results. 

Voltage profiles obtained with this sample are shown in figure 5 for selected 
values of applied voltage. The voltage probe sweeps across the sample about 
3 mm from the cover slide surface. 

When the segment is biased positively, the exposed interconnects collect elec- 
trons from the environment. When the applied bias voltage is less than +190 
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volts, the surface voltage profile shows that the bias voltage appears only at the 
interconnect regions with the cover slide surfaces remaining essentially at zero 
voltage. For thi applied voltage regime, the electron coupling current collected 
is proportional to the applied voltage and the area of the interconnects. 

When the applied voltage is greater than 190 volts, there is a transition in the 
collection phenomenon. There is a "snap-over" or sudden change in the 
surface voltage profiles. When this snap-over occurs, the bias voltage is no 
longer constrained to the interconnect region, but expands to encompass the 
cover slides as well. This results in essentially a Hat surface voltage across 
the array segment. The magnitude of this measured surface voltage is about 
50 volts less than the applied potential. The electron coupling current collected 
in this voltage regime is proportional to the applied voltage (less than 50 volts) 
and the array area. It is believed that this snap- over condition occurs due to 
secondary emission from the cover slides. The details of the theory to explain 
this phenomenon are being developed. 

When the array segment is biased negatively, a different phenomenon occurs. 

The voltage profiles indicate that ion collection is limited to the exposed con- 
ductor areas and that ion attractive potentials are significantly less than applied 
potentials. There is an additional phenomenon that occurs when the segment is 
biased negatively to a value greater than 700 volts-, the segment arcs (see fig. 1). 
This occurs on both the front and rear sides of the array (at times simultaneous 
front and back discharges were observed). This arcing occurred continuously 
during the testing of the segment at a rate of about 5 per minute. Testing at 
different plasma densities indicates that arcing occurs at different values of 
applied voltage, being higher at lower densities. At simulated geosynchronous 
altitude conditions this breakdown is expected at ~-8 to -10 KV based on pre- 
vious LeRC experiments. 


Experimental results presented here are preliminary. Trends observed have 

II 1 5 

been obtained with different samples in this ;uid other facilities ’ . Hence, 

it is believed that the array will collect electron currents as a function of ex- 
posed collector area and applied voltage up to positive values of 150 to 190 volts. 
Above this threshold, the current collection will be proportional to the panel 
area. At negative applied voltages, the ion current will lx* dependent only on 
the exposed conductor area. The detailed investigation of this phenomenon 
(both analytical and experimental) is continuing. 
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Simple Model tor Space System 


Calculations of array voltages and parasitic current collected as functions of 
array operating voltage were made for environments expected at 100 km and 
geosynchronous (both equatorial), using a simple model. The "space system" 
considered is that of figure a 20 meter x z50 meter solar array with a power 

o 

output of 500 kw ( -100 W m“). Some of the assumptions incoi*porated in the 
model are illustrated in figure 6. Two voltages for current collection are 
assumed, V and V . These are treated throughout this paper as magnitudes 
i.e. positive quantities. As figure 0(a) illustrates, V + + Y_ 1/2 V . Because 
the voltage gradient is assumed constant along the array, the length of the array 
available to collect ions or electrons is proportional to the ion or electron col- 
lecting voltage, i.e. f / f _ V + \ . V 4 and V_ are expected to be different 
because the electron and ion fluxes are different, as has been noted. 

The condition for the array to be in equilibrium with the surrounding plasma is 
that the net current to the array be zero, or. equivalently, that the magnitude of 
the electron current collected equal that of the ion current 

VV=jjA_ (1) 

Here, ami jj are the magnitudes of the electron and ion current densities 
striking the array, and A or A is the positive or negative array area col- 
lecting electrons or ions. 

The electron current density j is taken to be the random electron current den- 
sity. It is assumed that secondary electrons due to electron impact on the jxrsi- 
tive section of the array will be realtracted and thus contributed no net current 
density. 

The ion current density j. is taken to be the sum of two terms, one due to in- 
coming ions, and the other to secondary electrons from ion impact. The latter 
term is included because such secondary electrons will l>c repelled from the 
negative side of the array ami therefore add to the effective "ion" current. For 
purposes of this model, the current density due to these secondary electrons is 
taken to be 

l.tieJv - 

j. j. i » — 

,s Jl ° 1 4 10/V_ 

(V in kV) (2) 
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which gives a good fit to experimental data for protons normally incident on 

i r 

aluminum . Here, j Jo is the current density of incident ions and V is the 
negative array voltage, which corresponds to the energy of the ions on impact. 
While it is recognized that equation (2) is for protons rather than heavier ions 
(eg 0 4 ), it should give some estimate (within a factor 2 in the yield) of the effect 
of the secondary electrons on the current balance. 

Another consideration in the ion current density is the ram effect which occurs 
in low Earth orbit, where the velocity of the space vehicle through the plasma is 
greater than the thermal velocity of the ions. The situation for an orbiting solar 
array is sketched in figure 6(c). Assuming the array is oriented for normal in- 
cidence to sunlight, the velocity vector rotates around the array once per orbit. 
When 7 is normal to the plane of the array there is a ram ion current on either 
the cell side or the back of the array. In this case j. Q = jj ram . and ion col- 
lection takes place only on one side of the array. It is thus assumed that the ion 
density in the wake is sufficiently low that the contribution to the ion current 
from ions in the wake is negligible. When 7 is parallel to the array plane, 
only one edge sees the ram ions. As a limiting case of this condition, (assuming 
zero frontal area), j- o is Uiken to be just the thermal ion current for this 
orientation, and ion collection occurs on both sides of the array. That is, the 
collection is assumed to be from an isotropic plasma in this orientation. For 
either case, jj = j. Q + j.,_. 

There is no ram effect for electron collection, because the electron thermal 
velocities are much greater than the space system s oi*bital velocity. Thus, 
electron collection is always assumed to be from an isotropic plasma i. c. elec- 
tron collection occurs on both sides of the array. 


The areas for collection of electrons (A f ) and ions (A ) are taken to be 


A + = N + f + w eff+ 
A_ =N_f_w eff _ 


(•'0 


where the + or - subscript refers to the positive or negative voltage portions 

which collect electrons or ions. A is the area for particle collection, ( is the 

length of the panel for collection, w ^ is the effective width for collection and 

N is the number of sides of the panel collecting (N = 1 or 2). The effective 

17 

width w e £j is calculated from 

w off± = w + ' r s± < 4 > 
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where w is the panel width (20 meters) and r s is the sheath radius (see 
fig. 0(b)) which is calculated from the Child- Langmiur relation 



(electrons) 


(5a) 



(ions) 


(5b) 


and here, m t> or nij is the mass of the electron or ion, e the electronic charge 
(singly charged ions are assumed), e is the permittivity of free space (8. 85 x 
10 Farad/ meter), and V ^ the effective potential. 


V eff± = fV t 


(Hef. 17) 


( 6 ) 


where f is the fraction of the panel area occupied by exposed metal. In the 
present calculation equation (6) was used to define V^_ for all V_ since the 
experimental data indicate that the insulating surface in fact "shields" the 
voltage of the underlying cells for negative bias. Note that the data are for 
voltages I V_ I 1 kV. It is entirely possible that a voltage snapover could occur 
for larger negative voltages, or in the presence of lower energy ions. For 
collection of electrons, equation (6) was used to define V eff+ for voltages 
V < 200 volts. For V + s 200 volts, 

V efff = V + - 50 volts < 7 > 

based on the experimentally observed snapover in the electron collection data. 
The value of f in equation (0) was taken as 0.05 (i.e. 5 % exposed metal) for 
both voltage polarities (V, < 200 V), and for Ixith sides of the array. 


Equation (1) becomes 



with 
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isotropic collection 2 


j 


CO 


j, 
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f 0.05 V, 

LV 50 


ram collection 

V < 200 volts 

+ ~ 

V , < 200 volts 


w = 20 m 



f V 




f f + f _ = 250 m 


V + V 

i 



V 

op 


J 



(9) 


In the definitions of i and j. , T and T. are the temperatures of the 
J eo J io* e l * 

electrons and ions, respectively, and v is the velocity of the space system 
through the plasma. 


Equation (8) can be solved subject to conditions (9) to yield V + , V_, f, and f_. 
The electron and ion currents collected can then be determined. The parasitic 
current due to the plasma is at equilibrium, 

‘PLASMA = Je A + = ^i A - (10) 

This calculation was carried out for an environment characteristic of a 100 km 

l v 

orbit (see table I). Since v gti is greater than the ion thermal velocity, the 
calculation was done for both isotropic and ram collection conditions. Results 
are shown in figure 7. The line labeled 1 \H1LA A' indicates the current which 
is developed by a 500 kw array at operating voltage V . The parasitic cur- 
rent re P ros cnts a loss of current (and therefore power) available to 

the load. The curves indicate that about a 10% loss of current is expected for an 
operating voltage of 10 kV in the ram collection configuration, and a total loss 
of current by V Q p = 60 kV. For the isotropic collection condition, 10% loss 
occurs around V „ = 45 kV. A sun-oriented array will cycle between the iso- 
tropic and ram conditions twice per orbit (ignoring eclipse for the moment). 
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The apparent negative voltage shifts alxuit 10 r ,’ through the cycle, while the 
parasitic current shifts about an order of magnitude. 


The bars on the current curves in figure 7 indicate the region of V where 
the predicted negative voltage (V ) is in the range -750 to -1000 volts; this is 
the range of applied voltages which has been observed to result in arcing during 
laboratory experiments. According to the model, tins arcing range is reached 
well before the parasitic current drain becomes prohibitive. This arcing 
phenomenon is believed to be a technology problem in array design, and not to 
represent any absolute limit on array operating voltage. 


Those results indicate that the interactions of thermal plasmas with space sys- 
tems must be considered in the design of the large high voltage systems of the 
future. Both the parasitic currents and the arcing must l>e understood and tech- 
niques developed for eliminating their effects on space systems. As figure 7 in- 
dicates, it is not surprising that these effects have not presented spacecraft power 
system difficulties for systems flown to date. The highest operating voltage 
which has been used for spacecraft power systems is 100 volts, on Sky lab. The 
model predicts parasitic currents of 0.00 r (ram) and 0. 00:V' t (isotropic) for a 
100 volt system. Losses of this magnitude arc negligible. 


The calculation of floating voltages and parasitic currents was ;dso done for the 
quiescent geosynchronous environment (i.e. no substorm activity). The en- 
vironment parameters* 1 used are given in table II. Results of the calculation 
are shown in figure 8, with the results of the 100 km calculation drawn in for 
comparison. There is no ram effect in this case, since the ion thermal velocity 
is greater th;m the spac*cci*aft orbital velocity. The parasitic current at geo- 
synehronous is small (<5 r , of array current) for all operating voltages considered. 
There is the possibility of arcing for operating voltages >16 kV, based on the 
experimentally observed breakdowns at negative voltages of 8-10 kV for simu- 
lated geosynchronous conditions. 


There are several assumptions in the model used here which are worthy of note 
First, it lias been assumed that the array operating voltage is unaffected by any 
parasitic currents collected. In the real situation, the parasitic currents will 
shift the operating points of the solar cells, and cause the array voltage to fall 
until some equilibrium is reached. The amount of shift in operating voltage of 
a given string of solar cells will depend on its "position" in the array ("posi- 
tion" here refers to both physical location and electrical configuration) and the 
voltage- current characteristics of the cells. Because this loading effect has 
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not been included in the present calculations, the results should not be taken as 
a quantitative description of the parasitic currents collected; they merely give 
an indication of the range of design operating voltages which can be expected to 
be severely impacted by this class of interactions. 

Secondly, the modeling has assumed planar sheath growth (by using the Child 
Langmuir formula for r s ) throughout. This assumption becomes less valid as 
plasma density decreases and as voltage increases. Thus the model predictions 
are more reliable for the low Earth orbit case than for geosynchronous. In 
particular, the predictions at large operating voltages in geosynchronous orbit 
should be interpreted as order of magnitude estimates. 


Thirdly, no attempt has been made here to account for the effects of magnetic 
fields on the current collection. There are at least two sources of magnetic 
fields, the Earth's field and any space system fields, hi general, the presence 
of a magnetic field is expected to reduce the current that will be collected. It 
may be possible to "shield" the system by an appropriate configuration driving 
currents to create magnetic fields. This possibility should lie assessed during 
total system design and optimization. 


Finally, effects due to photoelectron emission and/or collection have not been 

included in the model. Its impact on the results cam, however, lx* estimated. 

-5 . o 

'Hie current density clue to photoemission should be of order 10 amps/m . 

This is small compared to electron and ram ion current densities at 100 km, 

and about equal to the random ion current. The maximum total photocurrent 

. - 5 

available for collection can be estimated as (f x w)xl0 0.05 amp. Since 
this is small compared to the plasma currents (except for the low voltage iso- 
tropic case) at -100 km, the calculation should not be affected greatly. In the 
case for the geosynchronous environment, the photocurrent density is greater 
than the thermal current density for both electrons and ions. It is therefore 
possible that the 1 photocurrent will dominate the equilibration, i.e. that the 
positive and negative voltages at equilibrium will be approximately equal. Since 
the total photocurrent is a 0.05 amp, the current collected may increase to that 
level; however, it still does not provide enough parasitic current to impact 
array operation seriously. The effect of shifting the relative magnitudes of the 
positive and negative voltages so that they were approximately equal would be 
to shift the * 1 expected onset of arcing" point towards higher operating voltages. 


The conelusion reached is, then, that the interactions of high voltage systems 
with the thermal plasma environment must be considered in the design of high 
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voltage systems for space applications. The possible detrimental effects to be 
avoided are excessive parasitic currents and arcing. While the calculations 
here have all been related to a hypothetical solar array, any high voltage com- 
ponent exposed to the environment can be expected to interact in the same gen- 
eral fashion. 

SUMMARY AND CONCLUSIONS 

Interactions of high voltage systems with the thermal plasma environment can 
give rise to parasitic current collection by high voltage systems (effectively a 
power loss mechanism) and to arcing for portions of the system sufficiently 
negative with respect to plasma potential. The simple analysis presented here 
indicates that the parasitic currents are of significant concern for systems with 
operating voltages of ~10 KV or greater in low Earth orbit. Parasitic currents 
are not expected to be of great concern for systems at geosynchronous (at least 
for operating voltages <100 KV, the highest considered here). The arcing of 
negatively biased segments of the system appears to be a potential problem for 
floating systems with operating voltages £l. 5 KV in low Earth orbit and with 
voltages >16 KV in geosynchronous. Design guidelines for high voltage systems 
are thus required. A program to investigate these phenomena and to develop the 
required design criteria is now underway at the Lewis Research Center. 

The interactions discussed here arc not believed to present insurmountable ob- 
stacles to future space systems. It is necessary however to design the new 
systems for space applications with an under standing of the ways in which their 
environments can affect them in mind. 
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Table 1 

400 KM ENVIRONMENT 

ions: O^g: ?j Q = 2xl0 11 /m* 
kT, = 0.03 eV 

l 

kT. 

= 0. 29 km/sec 

27rm. 

_ 4/3 

electrons: e : tj q = 2x10 /m 


kT ( , = 0.22 eV 


Spacecraft orbital velocity: 

7.7 km/sec 

Current densities: 


j. = 9. 4x10"^ amp/ m" 
J io 

(isotropic) 

jj Q = 2. 6x10” * amp/ m^ 

(ram) 

j„ = 2. 4x10' J amp/m 2 
J eo 

(isotropic) 


Table 2 

GEOSYNCHRONOUS THERMAL ENVIRONMENT 
ions: H : tj = 9. OxioVm’* 
k r L = 1.35 eV 
kT. 

= 4. 53 km/ sec 

27rmj 

- 7 3 

electrons e~: rj = 9. 5x10 /m 

kT = 1.35 eV 

Spacecraft orbital velocity: 3 km/sec 

Current densities: 

-8 2 

j. Q = 7x10 ' amp/m (isotropic) 

0 2 

j = 3x10 ° amp/ni (isotropic) 


! 



M5 ». 



’? V<jOR 
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Figure 2. - Orbital charged particle environments. 
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Figure 3. - Large solar array interaction model 500 kW array. 
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Figure 4. - Pictures o( array (composite*. 
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Figure 5. - Experimental results. 
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Figure 6. - Interaction computation models. 



Figure 7. - Plasma/solar array coupling currents vs solar array oper- 
ating voltage for a 500 kW array-, in thermal environments: 400 km 
environment. 
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Figure 8. - Plasma/solar array coupling currents vs solar array oper- 
ating voltage for a 500 kW array; in thermal environments: 400 km 
vs geosynchronous. 
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